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ABSTRACT 


Naval Supply Systems Command (NAVSUP) supports Navy, Marine Corps, Joint 
and Allied Forees with their inventory of more than 430,000 items worth $21 billion 
using several distribution sites. Choosing the optimum order-point and order-quantity 
for eaeh item is important to meet the stoehastie demand while satisfying multiple 
restrietions sueh as budget and maximum number of orders. The Site Demand-Based 
Level Inventory Optimization Model (SIOM) is a mixed-integer, linear program 
developed at the Naval Postgraduate Sehool to provide NAVSUP planners with guidanee 
on this eomplex problem. Ongoing tests have been sueeessful, but SIOM’s eomputational 
run times are long. This thesis introduees a new, faster reformulation (SIOMsQ) that 
approximates the solution of the same problem by redueing the possible eandidate sets of 
order-points and order-quantities for each item. We find that the solutions suggested by 
SIOMsQ are better than or very close to those of SIOM in test cases provided by 
NAVSUP, with substantially shorter computational times. Therefore, we recommend 
using SIOMsQ versus SIOM. 
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EXECUTIVE SUMMARY 


The mission of the Naval Supply Systems Command (NAVSUP) is to provide the 
Navy and Joint Warfighter with global logisties and quality-of-life support (U.S. Naval 
Aeademy 2016). NAVSUP Weapon Systems Support (NAVSUP WSS) is responsible 
for a large fraetion of NAVSUP’s National Item Identifieation Numbers (NIINs) with 
more than 430,000 items (NAVSUP 2016). As its primary mission, NAVSUP WSS has 
to provide all the Navy and Marine Corps units with the items they need with minimal 
or no delay. 

In order to maximize the operational readiness of hundreds of units, NAVSUP 
WSS needs to keep adequate wholesale and retail inventory levels for all items subjeet to 
multiple eonstraints. As a measure of effeetiveness, “fill rate” is used to speeify the 
expeeted fraetion of orders for whieh a replaeement is available (i.e., on hand when 
demand oeeurs). Aehieving a good fdl rate requires ehoosing an order-point, 5, and an 
order-quantity, Q, for eaeh item. 

In order to find optimum order-points and order-quantities, Salmeron and Craparo 
(2016) have developed Site Demand-Based Level Inventory Optimization Model (SIOM) 
to guide stoek level deeisions for NAVSUP. Beeause of the large seale and non-linear 
aspect of the problem, SIOM pre-generates candidate order-quantities for every NUN, 
and then optimizes order-points and order-quantities (restricted to pre-generated ones). At 
the time of this research, SIOM is undergoing testing for ten test cases at NAVSUP WSS. 
Preliminary results are promising, albeit computational times are long. This thesis 
develops SIOMsQ, an efficient reformulation of SIOM which improves both solution 
quality and computational time. 

SIOMsQ approximates the solution of SIOM by reducing the possible candidate 
sets of order-points and order-quantities for each item. That is, we proceed as in SIOM 
but limit the values of, not only Q, but also s, to a list of pre-generated candidate {s, Q) 
pairs. For example, SIOMsQ can generate up to 10 candidate order-quantities for every 
NUN, and then 20 candidate order-points for each order-quantity and NUN. First, we set 
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upper and lower limits of both 5 and Q using the data provided by NAVSUP (shelf life, 
maximum stock level, average monthly demand, etc.). Then we develop a strategy to 
select the candidate order-points and order-quantities between their respective limits, 
such that the differences between consecutive candidates are approximately the same. 

SIOMsQ is formulated as a pure integer, linear program that chooses exactly one 
candidate for each item, so that the overall choice is feasible for all items. The objective 
of SIOMsQ is to minimize the total expected penalty by choosing the optimum (already 
generated) {s, Q) pair for each NIIN while satisfying the constraints such as budget, 
maximum number of orders per month and shelf life among others. Since, for each 
generated {s, Q) pair, we can pre-calculate the associated fdl rate, cost, penalties, and 
other data that would apply should such pair become the decision adopted for the item, 
the resulting model does not need the complex constructs of SIOM. The possibly large 
number of binary decision variables due to candidate pairs is compensated by a reduction 
in model complexity, as proven by our computational results. 

SIOMsQ is run in General Algebraic Modeling System (GAMS) using CPLEX 
solver (GAMS 2016). We run SIOMsQ with various candidate (5, Q) pair sizes which 
range from 7 to 20 candidate Q-values, and similarly for candidate 5-values for each Q. 

For all of the ten test cases (from real-world problems provided by NAVSUP), the 
solutions SIOMsQ achieves are better than or very close to those of SIOM, but the 
running times of SIOMsQ are substantially shorter. The improvements are most 
compelling in large-scale problems, where SIOM has notable solving difficulties. 

As a conclusion of this research, we recommend that running SIOMsQ using 
either 7 Q-values and 7 5-values for each item (if very fast runs with acceptable solution 
quality are required), or 10 Q-values and 20 5-values for each item (if solution quality is 
very important, still with acceptable running times). 
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I. 


INTRODUCTION 


A. BACKGROUND 

The mission of the Naval Supply Systems Command (NAVSUP) is to provide the 
Navy and Joint Warfighter with global logistics and quality-of-life support (U.S. Naval 
Academy 2016). Involving about 345 military personnel and more than 24,000 civilian 
employees, NAVSUP holds inventory for approximately $21 billion and spends a yearly 
budget of $3.5 billion (NAVSUP 2016a, b). 

The mission of NAVSUP Weapon Systems Support (NAVSUP WSS) is to 
deliver weapon systems supply support to the Navy, Marine Corps, Joint and Allied 
Forces (NAVSUP 2016c). NAVSUP WSS is responsible for a large fraction of 
NAVSUP’s National Item Identification Numbers (NIINs) with more than 430,000 items 
(NAVSUP 2016a). 

As its primary mission, NAVSUP WSS has to provide all the Navy and Marine 
Corps units with the items they need with minimal or no delay. 

In order to maximize the operational readiness of hundreds of units, NAVSUP 
WSS needs to keep adequate wholesale and retail inventory levels for all items subject to 
multiple constraints. A key complication in the analysis is the fact that demands for the 
NIINs are stochastic. 

As measure of effectiveness, “fill rate” is defined as the expected fraction of 
orders for which a replacement is available (i.e., on hand when demand occurs). A 
detailed formulation is included in Section II.C. 

Due to limitations such as inventory budget, number of orders per month, storage 
space and shelf life, it is not possible to achieve a perfect fill rate of 100% for all NIINs. 
Instead, NAVSUP WSS specifies minimum required fill rate levels for all items, based 
for example on the importance of the item to readiness. 

This thesis develops a mathematical optimization model to guide inventory 
decisions for NAVSUP. 
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B. REVIEW OF INVENTORY MODELS 


Silver et al. (1998) point out the importance of inventory management and 
production planning and scheduling of several industries and organizations in economic, 
medical and military aspects. When modeling the inventory management system of 
individual items with probabilistic demand, Silver et al. introduce two techniques: 
continuous review and periodic review. 

In continuous review, the state of the inventory is assumed to be known at all 
times; for example, an order can be placed at any point in time. In periodic review, the 
inventory is verified (and related decisions are made) at discrete points in time (e.g., on a 
given weekday). For this study, we are solely concerned with continuous review models. 

Silver et al. (1998) describes two types of continuous review inventory 
management systems: order-point, order-quantity {s, Q) system and order-point, order- 
up-to-level {s, S) system. (Note: order-point and order-quantity are sometimes referred to 
as reorder point and order size, respectively. We use the same terminology as in Silver et 
al. (1998) in the rest of this thesis.) 

A backorder is a demand that cannot be filled at the time it occurs (i.e., an 
outstanding order for which the customer will need to wait). Net stock is defined as the 
number of items on hand, whereas inventory position includes the number of items 
ordered and backordered in addition to net stock. Inventory position at any given point in 
time can be formulated as follows: 

Inventory Position = Net Stock + Orders - Backorders 

In an {s, S) system, when the inventory position drops below 5, an order is placed 
such that the inventory position reaches the order-up-to level, S. Different amounts of 
items can be ordered at different times depending on how far the inventory position is 
away from the order-up-to-level. On the other hand, in an (s, Q) system, when the 
inventory position drops below 5, a pre-determined order amount Q is placed. 

It should be kept in mind that an order is placed when the inventory position 
drops below order-point, 5. If orders were placed according to the net stock, there could 
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be unnecessary orders that stem from not waiting the lead time of the previous order 
(Silver et al. 1998). 

In Figure 1, inventory position is represented by dashed line and net stock or both 
the inventory position and net stock are shown by the solid line. There are demands of 
different amounts at several points in time. When the inventory drops below s at time A, 
an order is placed and the inventory position increases by Q. L represents the order lead 
time, which is assumed to be deterministic. 
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Figure 1. is,Q) Inventory System. Source; Silver et al. (1998). 

For a certain item, choosing a large order-point provides a good fill rate, but this 
increases the amount of money tied up in inventory. Also, some items might exhaust their 
shelf life before they are demanded. On the other hand, if a small order-point is selected, 
it could cause backorders, thus reducing the fill rate. 

Similar to the order-points, order-quantities have an important effect on fill rates. 
If a small order-quantity is selected for an item, backorders could arise depending upon 
the lead time (the time elapsed between an order placement and its arrival) of the item. 
Also, a small order-quantity requires orders to be placed more frequently, incurring 
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additional costs and preventing potential savings sueh as discounts from seller. On the 
other hand, a large order-quantity could cause some items to perish before demanded. 

Determining the order-points and order-quantities for thousands of items in order 
to achieve target fill rates is a complieated task because all the NIINs share a budget limit 
and an order limit (the number of order doeuments that can be placed in a given time 
period). 

For the problem at hand, an ( 5 , S) model is computationally intractable; therefore 
we develop an (5, Q) model and approximate the optimal (5, S) solution by setting 
S = s + Q. 

Salmeron and Craparo (2015) introduee the Wholesale Inventory Optimization 
Model (WIOM). Given the order-quantities, WIOM optimizes the order-points of all 
items simultaneously by minimizing the expected weighted deviation from the target fill 
rates of each item. 

Roth (2016) eompares the performances of three different inventory management 
models using a simulation: simple caleulation of fill rates, SPO and WIOM. His results 
show that the order-points suggested by WIOM provide higher fill rates than the two 
other models. 

The key idea this thesis exploits is the enumeration of possible values of deeision 
variables to enable pre-proeessing of information with whieh to build an approximate 
optimization model. Enumeration is a common technique in optimization. For example, 
Kolodziej et al. (2013) use this method in optimization of the multi-period blend 
scheduling problem. In order to overeome the diffieulties of the non-eonvex MINLP 
problem, they enumerate to approximate the optimum blending fractions of various 
products. 

C. PROBLEM STATEMENT 

OPNAVINST 4441.12D defines demand-based items as "items that have a 
relatively high issue rate" (Office of the Chief of Naval Operations 2012). 
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NAVSUP WSS is now using Morris Cohen Associates' ‘Service Planning and 
Optimization (SPO)’ tool for planning order-points and order-quantities for 
demand-based items. However, they do not document the mathematical model or 
algorithms used to determine such values (Roth 2016). 

In order to find optimum order-points and order-quantities, Salmeron and Craparo 
(2016) have developed Site Demand-Based Level Inventory Optimization Model 
(SIOM). Because of the large scale and non-linear aspect of the problem, SIOM 
pre-generates candidate order-quantities for every NUN, and then optimizes order-points, 
5, and order-quantities, Q (restricted to pre-generated ones). SIOM is implemented using 
the General Algebraic Modeling System (GAMS 2016a) with the GAMS CPLEX 
(GAMS 2016b) as solving engine. For completeness, the formulation of SIOM is given in 
the Appendix. 

At the time of this research, SIOM is undergoing testing at NAVSUP WSS. 
Preliminary results are promising, albeit computational times are long due to the 
complexity of the resulting optimization model. 

In this study, we introduce a new technique based on pre-generation of both 5 and 
Q candidates. That is, we proceed as in SIOM but limit the values of, not only Q, but also 
s, to a list of pre-generated candidate (s, Q) pairs. The new optimization model is referred 
to as SIOM with pre-generated s-Q (SIOMsQ). 

The motivation for enumerating (s, Q) pairs by pre-generation is two-fold; First, 
computational time of SIOMsQ is expected to be shorter due to using a restricted set of 
candidate {s, Q) pairs. Second, and possibly the most important, once {s, Q) pairs are 
generated, fill rates can be calculated (for each item and (s, Q) pair) using whichever 
method is deemed appropriate, such as closed-form equations, simulation, etc. For 
example, Roth (2016) has devised an accurate simulation of fill rates, which could 
potentially be used in SIOMsQ as input data. In SIOM, since the order-points are 
decision variables themselves, the possibilities to calculate fdl rates are reduced to an 
approximation via a series of equations and constraints embedded in the model, making it 
substantially more complex. For comparison purposes, we only calculate the fill rates for 


5 



the pre-generated ( 5 , Q) pairs, mimicking the values that would have been calculated if 
SIOM had been used. 

D. RESEARCH OBJECTIVES AND SCOPE 

The objective of this research is to develop an algorithm to generate reasonable 
(s, Q) pairs for every NUN, and to create the SIOMsQ formulation to approximate the 
solution of the original SIOM faster. SIOMsQ seeks to strike a proper balance between 
solution quality and solvability of the problem. That is, SIOMsQ may relinquish solution 
quality for a faster solution. 

This research assumes the existing SIOM release 1.2.2 is a valid approach to the 
site demand-based inventory problem. 

Some distribution sites have two types of NIINs and they are treated as different 
sites. NMC stands for “Navy-Managed Consumables” and BP28 refers to “Navy 
Working Capital Fund Budget Program 28” (Salmeron and Craparo 2015). The SIOMsQ 
approach will be compared to SIOM in the following test cases already provided by 
NAVSUP: 

1. Norfolk (NMC) 

2. Norfolk (BP28) 

3. Yokosuka (NMC) 

4. Yokosuka (BP28) 

5. Kings Bay (NMC) 

6. Kings Bay (BP28) 

7. San Diego (NMC) 

8. Bangor (BP28) 

9. Jacksonville (BP28) 

10. Key West (BP28) 
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II. METHODOLOGY 


A. SIOM CHARACTERISTICS 

SIOM (Salmeron and Craparo 2015) minimizes the aggregate nonlinear penalty 
for not satisfying the required fill rates of each item. The model is non-separable by item 
due to constraints on total budget and maximum number of orders per month, each of 
which involves all NlINs in the problem. 

Due to non-linear relationships between order quantities, fill rates and number of 
orders placed, SIOM pre-generates a list of candidate Q-values for every NUN, which 
allows us to formulate SIOM as a mixed-integer linear optimization problem. 

The order quantities (currently up to ten candidate Q-values for every NUN) are 
generated using the following ideas: 

• Set a reasonable lower and upper limit for Q, based on the data (average 
demand, lead time, shelf life, etc.) already provided by NAVSUP WSS, 
and 

• Select Q = l as a candidate value and then remaining order-quantities 
between the upper and lower bounds (inclusive) such that step sizes 
between two consecutive candidates are approximately the same. 

After generation of order-quantity candidates, SIOM optimizes s and Q (as 
restricted) for all NllNs. 

However, in some of the sites that have thousands of items, SlOM’s difference 
between the upper and lower bounds of the solution are large even after 5 hours of 
computation. 

To be able to obtain a faster solution, SIOM incorporates a grouping heuristic 
algorithm (GHA). In GHA, items in a certain distribution site are separated into smaller 
groups, constraints are adjusted proportionally and sub-problems are solved individually. 
At the end, objective values of sub-problems are summed and the approximate solution is 
achieved. For intermediate size sites (3,000-5,000 items), the grouping algorithm gives 
good results in a reasonable amount of time. In larger cases (5,000+ items), in order to 
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obtain an acceptable solution, the number of groups should be inereased whieh 
potentially reduees the quality of the solution. 

B. SIOMSQ DEVELOPMENT 

1. Definitions 

The development of SIOMsQ needs the following terms to be defined: 

Shelf life ): Shelf life for a given item (in months). 

Maximum Stoek Level (MSL): Inventory up-to-level for a given item (i.e., same 
as S. = s,+Q^). 

Expeeted monthly demand {xf ): Expeeted monthly demand for a given item. 

Allowanee (^ ): Minimum MSE level required for a given item. If not speeified, 
the default value is 0 for eaeh item. 

Maximum Months of Supply (S ): A parameter that restriets the MSE of eaeh 
item. Eor instanee, if S is 12, MSE of an item eannot exeeed 12 times the expeeted 
monthly demand of the item. 

Minimum Months of Supply (S): The parameter that sets a lower bound on MSE 
of eaeh item. If S is (for example 0.5), maximum stoek level eannot be less than 0.5 
times the expeeted monthly demand of the item. 

Minimum order-point (^.): Eower bound on 5 for a given item. Default is -1, 

whieh means that we would wait until there are no items on-hand and a new demand 
arrives to plaee an order of the item. 

2, Generation of Parameters s and Q 

Order-points and order-quantities are integer deeision variables whieh ean take 
potentially many values. To simplify the problem, SIOMsQ uses a redueed set of 
eandidates for both order-points and order-quantities. Eor example, in order to generate 
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candidate ( 5 , Q) pairs we ean generate up to 10 eandidate order-quantities for every NUN, 
and then 10 candidate order-points for eaeh order-quantity and NUN. 


a. Generation of Order-Quantities 

To generate order-quantities, we preserve the same seheme as in SIOM (Salmeron 
Craparo 2016); 

• “1” is always the first candidate Q for each item. 

• Set lower bound on the remaining Q values ( Q.) as max {S x , 2} . 

• If shelf life of the item is less than S , set upper bound on Q {Q .) as 

X xf . Otherwise, set Q. as max{^^., S xxf}. 

• Set Q. and as the seeond and the last eandidates, respectively. 

• Seleet the remaining eandidates between Q. and Q ., sueh that the 

differenees between eonseeutive eandidates are approximately the same 
(i.e., round to the nearest integer for fraetional numbers). 

b. Generation of Order-Points 

In SlOMsQ, generation of eaeh eandidate 5 depends on already generated 
eandidate order-quantities. Therefore, every eandidate Q has its own set of order-points: 

• “-1” and “0” are always the first two eandidate order-points for each item. 
Candidate “-1” stands for baekorders. 

• Set lower bound on the remaining 5 values (y,) as max{^. - Q ., 1} . 

• Set upper bound on 5 ( 5 ,.) as min{5,.^ x xf -Q^,S x xf + Q;}. 

• Set ^ and s^ as the third and the last eandidates, respeetively. 

• Seleet the remaining eandidates between ^ and 5 ,, sueh that the 

differenees between consecutive candidates are approximately the same 
(i.e., round to the nearest integer for fraetional numbers). 

C. MATHEMATICAL MODEL 

1. Indices and Index Sets 


i 


Item (i.e., NUN), for i g I. 


9 



h Index of candidate ( 5 , Q) pairs for a given item i, /z e {1,2,...}. 

m Index for the penalty segments of fill rates (e.g., m^M = {1,2,5}). 

/ Group of items, for I e L. 

2. Input Data and Parameters [Units] 

f, Lead time for item i. [quarters] 

X, Expected demand for item / during the lead time, [units of issue/lead 
time] 

h. Number of candidate pairs (order-point and order-quantity) for item /. 

s-i^ 5-value (order-point) of the h-th candidate for item i. [units of issue] 

2 -value (order-quantity) of the h-th candidate for item /. [units of issue] 

^ Lower bound on order-point of item i (default is -1). [units of issue] 

Upper bound on order-point of item /. [units of issue] 

Allowance for item /. [units of issue] 

A- 1 if allowances are activated, and 0 otherwise. 
c. Cost per unit of item /. [$/unit] 

r Maximum number of total expected orders per month for all items, 
[orders/month] 

Si^ Shelf life for item /. [quarters] 

S Maximum months of supply for any item, [months] 

Initial order-point used to enforce persistence for item /. [units of issue] 

S[' Penalty for deviation from initial order-point for item /. [unitless] 
y'’ Overall persistence weight, [unitless] 
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S,^ Penalty for deviation of item i above maximum months of supply, 
[unitless] 

Overall weight for MSL deviations above maximum months of supply, 
[unitless] 

1. Group of item /. [unitless] 

fj Required fill rate for any item in group I. [fraetion, unitless] 

w, Weight for aehieving the required fill rate of any item in group I. 
[weight units] 

3. Decision Variables 

r, ,, 1 if eandidate pair h is selected for item /, and 0 otherwise. 

Q Order-quantity for item i. [units of issue] 

Si Order-point for item /. [units of issue] 

4. Derived Data 

Fill rate of item i for candidate pair h . [fraction, unitless] 

Penalty for selecting candidate pair h for item i. [unitless] 

s, ]; Deviation down with respect to initial order-point for item i and candidate 

pair h . [units of issue] 

s' Deviation up with respect to initial order-point for item i and candidate 
pair h . [units of issue] 

S*,^ Deviation for MSL of item i and candidate h above the average demand 
during the maximum months of supply for the item, [months] 
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c-f^ Expected number of cycles per lead time period for item i if candidate 

pair h is selected, [unitless] Note; A cycle is defined as the time period 
between consecutive orders, or the lead time; whichever is shorter. 

B.f^ Expected fraction of backorders during a cycle for item i if candidate 
pair h is selected, [fraction, unitless] Note: The calculation method 
appears in Equation (1). 

X'.f^ Random variable for the demand of item i in a cycle if candidate 
pair h is selected, [units of issue]. 

Deviation from the required fill rate for item i if candidate 
pair h is selected, [fraction, unitless] 

Deviation penalty for item i , in fill rate penalty segment m if candidate 
pair h is selected, [unitless] 

Maximum deviation allowed for fill rate of item i , in penalty segment m. 

2 

- ffl 

(Calculated as = /, —7 ) [fraction, unitless] 

' 

j^M 

^ Penalty for deviation from the required fill rate for items of group I in 
penalty segment m. (calculated as w, ^ = mwj ) [unitless] 

After generating all the candidate ( 5 , Q) pairs for all items, we can pre-calculate 
the fill rate of each item for each candidate pair, h = (s, Q). 

The only reason that could cause fill rate to be less than 100% is backorders. 
Therefore, in order to calculate the fill rate of every item for every corresponding h, first 
we need to compute the expected fraction of backorders, , during a cycle (see 
definition of ) for each item and for every h. The expected fraction of backorders can 
be calculated using the following formula: 
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1 

4=^j(x-4)fx,^(x)dx 

^ih S:U 


( 1 ) 


In Equation (1), stands for the random variable for the demand in a eyele with 
probability distribution given by density funetion (or mass function) f^^, (x) and 


X. 


s'ih = ^ih - (^ih - 1) Qih ^ where c,, = max 11, ^. 


After calculation of expected fraction of backorders, fill rate of each NUN can be 
computed (for every h) as follows: fih=l- B.,^ . This formula can return f.^ < 0 in which 
case we use /* = 0 . 


Penalties are applied for the items that cannot achieve the required fill rate. For 
every site, there are certain groups of items that all have the same required fill rates. 
Deviation from the required fdl rate is formulated as - max{0, - fn,}, where fi is 
the required fill rate for item i of group I. 

Before calculating the total penalties for fill rate deviations, we need to compute 
the deviation penalties in every penalty segment for each item, , which is calculated 
as follows: 


f w 

J im IjM ’ 


if ft > I 




m'<m 


f _ ^ 




r dev X ’ r - 

J ih J im' 


if Z ff, < ft 

^ y /■ 

J im 




m'<m 

0, 


if/;""i Z 



m<m-\ 


Finally, the deviation penalty for selecting candidate pair h of item i can be 
calculated as = w, . 


Another derivation that has to be computed as an input to SIOMsQ is s* and . 
They are calculated as 4 = max{s“-s,^,0} and s“ = max{s,,,-s",0}. 
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Lastly, in order to calculate MSL deviation penalty, s* has to be derived as 
S* = max(s,., + Q,„) - S,0}. 

X. 


5. SIOMsQ Formulation 

SIOMsQ can be stated as the following mixed-integer, linear problem: 
SIOMsQ: minX 

’ ’ i h\h<h: ?■ h\h<h: n*l.J j li\h<h: “i” J- 


( 2 ) 


Subject to: 


Qi=Yj_Qih^ih Vi 

h\h<hi 

(3) 

h\h<hi 

(4) 

h\h<h, 

(5) 


(6) 

i h\h<h^ 

(V) 

+ V/ 

(8) 

V/ 

(9) 

g-> 0 and integer v/ 

(10) 

i'; > -1 and integer V/ 

(11) 

r,, £{0,1} yi,h\h<h 

(12) 


The goals of the objective function (2) are to minimize (a) deviations from target 
fill rates for all items (larger penalty rates for being far away from target fill rate); (b) the 
penalties for deviating from current order-points; and (c) the penalties for exceeding 
maximum months of supply. Persistence penalties and MSL deviation penalties are 
optional and can be removed from the objective function by setting y'’ = 0 and = 0 

9 

respectively. 

Equations (3) - (5) restrict the model in choosing only one ( 5 , Q) pair for an item 
and set the values of selected 5 and Q. 
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Equation (6) establishes a budget restrietion on the MSL eost of ah items. 
Equation (7) restriets the expeeted number of orders plaeed in a month from exeeeding a 
given limit. 

Equation (8) limits the lower and upper bounds of MSE. Stoek levels should be 
greater than the ahowanee and less than the minimum amount that would eause items to 
perish (under expeeted demand assumptions). 

Equations (9) - (12) set the domains of the deeision variables and parameters. 

6. Two-Step Approach 

While running SIOMsQ (regular method), we empirieahy found that it takes 
substantial time for GAMS/CPEEX to identify an initial solution before starting 
braneh-and-bound. In faet, in rare oeeasions, a feasible solution (i.e., integer) eannot be 
found in the allotted time of two hours. In order to speed up SIOMsQ solution, we have 
developed the following two-step approaeh; 

• In the first step, we solve a further restrieted SIOMsQ, eahed SIOMsQ'. 
This uses only a random subset from the original eandidate (s, Q) pairs, 
for eaeh NUN. SIOMsQ' shortens the time for an initial solution 
ealeulation. We do not fully solve SIOMsQ': we stop after the first 
feasible solution has been identified. We denote by h'. < h ., the number of 
eandidate {s, Q) pairs for item / used in SIOMsQ'. 

• In the seeond step, we run the full SIOMsQ (with ah eandidate pairs) 
using the initial solution found by SIOMsQ'. 

D. COMPUTATIONAL EXPERIENCE 

Both SIOM and SIOMsQ are exeeuted on a Eenovo W520 laptop eomputer with 
16 GB of RAM and two 2.5 GHz Intel Core i7 proeessors. We direet GAMS/CPEEX to 
use up to four parallel threads, whieh speeds up ealeulations by a faetor of 2.0 in some 
eases. Ah SIOMsQ eases are solved to within 1% of optimality or using two hours of 
eomputational time, whiehever eomes first. SIOMsQ' is solved until the first feasible 
solution is identified. 
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III. COMPUTATIONAL ANALYSIS 


We run SIOMsQ for ten eases provided by NAVSUP. We present the results of 
these runs with various options and eompare them to SIOM. 

A. RUNNING OPTIONS FOR SIOM AND SIOMSQ 

SIOM is executed with time limit of 2 hours or 5 hours (if 2 hours is not enough 
to decrease the optimality relative tolerance to 1%) for each site. SIOMsQ always has up 
to 2 hours of running time. For the two-step approach, the time spent in SIOMsQ' is not 
counted in the two-hour limit allowed for the full SIOMsQ. Additional time is allowed 
for pre-processing of the data. 

Apart from full group size, SIOM may incorporate a grouping heuristic in which 
items are grouped into either 500 or 1,000 NIINs in order to solve the problem faster. The 
grouping algorithm is used only for the cases where an acceptable solution is not 
achieved within given time limit. SIOMsQ is always run with all NIINs. 

Both SIOM and SIOMsQ use 1% optimality relative tolerance. That is, while 
branch-and-bound is being executed, when the relative difference between the incumbent 
objective function value and the best possible value of the solution (lower bound) drops 
below 1%, the solver stops and the incumbent solution is reported. 

For a given set of pre-generated Q-values, lower bound on the objective function 
value is only guaranteed to be valid for SIOM if run with no grouping of NIINs. Lower 
bounds provided by SIOM with a grouping heuristic and by SIOMsQ are only estimates. 

SIOM pre-generates #2 = 10 candidate order-quantities and fully optimizes 
order-points for each item. In SIOMsQ, however, both s and Q candidates are 
pre-generated. In order to see the effects of number of pre-generated order-points and 
order-quantities on the solution quality and running time, we use various sets of {s, Q) 
pairs. We denote by {#s, #Q) a specific run where we pre-generate #Q candidates for Q 
for each NUN, and then #s candidates for s for each of the Q candidates. We test the 
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following combinations of (#s, #Q): (7, 7), (7, 10), (10, 10), (15, 10), (15, 15) and (20, 

10 ). 

SIOMsQ is run with the regular method (in one step without pre-generating a 
feasible initial solution), and with two-step approaeh (deseribed in Seetion II.C). In the 
latter ease, we use SIOMsQ' with a redueed number of eandidate pairs h,' = 50 V / in 
order to generate an initial solution faster. 

B. RESULTS 

1. Yokosuka (NMC) 

Yokosuka (NMC) is the smallest ease with 35 NIINs. SIOM solves the problem in 
less than a minute with no penalty applied to the objeetive funetion value. That is, the 
required fill rates are satisfied optimally with the solution suggested by SIOM for eaeh 
item. (See Table 1.) 


Table 1. Yokosuka (NMC) SIOM Results 


Time Limit 

Group Size 

Objeetive 

Lower Bound 

Run Time (mins) 

2 hours 

fiill 

0 

0 

<1 


SIOMsQ also solves the problem in less than a minute for all possible solving 
options. However, when the number of 5 and Q eandidates are (7, 7) and (7, 10) the 
objeetive funetion value is slightly sub-optimal. The regular method and the two-step 
approaeh do not make any differenee in the solution for this ease. (See Table 2.) 

For this ease, the SIOMsQ optimization problem size ranges from 5,217 integer 
variables and 248 single equations in the (7, 7) ease to 5,708 integer variables and 248 
single equations in the (15, 15) ease. 
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Table 2. Yokosuka (NMC) SIOMsQ Results 


# 5 

#Q 

Method 

Objective 

Estimated 
Lower Bound 

Run Time (mins) 

7 

7 

Regular 

0.008 

0.008 

<1 

Two-step 

0.008 

0.008 

<1 

7 

10 

Regular 

0.008 

0.008 

<1 

Two-step 

0.008 

0.008 

<1 

10 

10 

Regular 

0 

0 

<1 

Two-step 

0 

0 

<1 

15 

10 

Regular 

0 

0 

<1 

Two-step 

0 

0 

<1 

15 

15 

Regular 

0 

0 

<1 

Two-step 

0 

0 

<1 

20 

10 

Regular 

0 

0 

<1 

Two-step 

0 

0 

<1 


2. Key West (BP28) 

Key West (BP28) is one of the smallest test eases with 53 items. SIOM gives a 
solution of 1.40 in less than a minute. (See Table 3.) 


Table 3. Key West (BP28) SIOM Results 


Time Limit 

Group Size 

Objective 

Lower Bound 

Run Time (mins) 

2 hours 

Ml 

1.40 

1.39 

<1 


The running time of SIOMsQ for Key West (BP28) site is also less than a minute. 
When it is run with 7 candidates of s and Q each, the solution is worse than all the other 
cases. This is an expected outcome, since coverage in 5 and Q is low. 

In the case of 15 candidates for both s and Q, the solution is even better than the 
lower bound suggested by SIOM. The reason is that SIOM uses 10 pre-generated 
order-quantities. Due to a better coverage in Q, the solution quality can potentially be 
better in the (15, 15) case. (See Table 4.) 
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Table 4. Key West (BP28) SIOMsQ Results 


# 5 

#Q 

Method 

Objeetive 

Estimated 
Lower Bound 

Run Time (mins) 

7 

1 

Regular 

1.47 

1.47 

<1 

Two-step 

1.47 

1.47 

<1 

7 

10 

Regular 

1.40 

1.40 

<1 

Two-step 

1.41 

1.40 

<1 

10 

10 

Regular 

1.40 

1.40 

<1 

Two-step 

1.40 

1.40 

<1 

15 

10 

Regular 

1.40 

1.40 

<1 

Two-step 

1.40 

1.40 

<1 

15 

15 

Regular 

1.38 

1.37 

<1 

Two-step 

1.38 

1.37 

<1 

20 

10 

Regular 

1.40 

1.40 

<1 

Two-step 

1.40 

1.40 

<1 


3. Kings Bay (NMC) 

There are 251 items in Kings Bay (NMC). SIOM aehieves 30.7 as the objeetive 
funetion value in less than a minute with the full group size. (See Table 5.) 


Table 5. Kings Bay (NMC) SIOM Results 


Time Limit 

Group Size 

Objeetive 

Lower Bound 

Run Time (mins) 

2 hours 

Ml 

30.7 

30.7 

<1 


When we run SIOMsQ with various sets of (5, Q) pairs, we always get the same 
objeetive funetion value of 30.7. Using regular method or two-step approaeh does not 
ehange the solution. (See Table 6.) 
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Table 6. Kings Bay (NMC) SIOMsQ Results 


# 5 

#Q 

Method 

Objective 

Estimated 

Lower Bound 

Run Time (mins) 

7 

7 

Regular 

30.7 

30.7 

<1 

Two-step 

30.7 

30.7 

<1 

7 

10 

Regular 

30.7 

30.7 

<1 

Two-step 

30.7 

30.7 

<1 

10 

10 

Regular 

30.7 

30.7 

<1 

Two-step 

30.7 

30.7 

<1 

15 

10 

Regular 

30.7 

30.7 

<1 

Two-step 

30.7 

30.7 

<1 

15 

15 

Regular 

30.7 

30.7 

<1 

Two-step 

30.7 

30.7 

<1 

20 

10 

Regular 

30.7 

30.7 

<1 

Two-step 

30.7 

30.7 

<1 


4. San Diego (NMC) 

San Diego (NMC) has 468 NIINs in its inventory. SIOM solves the problem in 19 
minutes and the objeetive function value is 103. (See Table 7.) 


Table 7. San Diego (NMC) SIOM Results 


Time Limit 

Group Size 

Objective 

Lower Bound 

Run Time (mins) 

2 hours 

Ml 

103 

102 

19 


We obtain approximately the same objective function value as the solution in 
SIOM when we run SIOMsQ for any number of candidate s and Q or solution method, 
but it takes less than a minute in all cases. (See Table 8.) 
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Table 8. San Diego (NMC) SIOMsQ Results 


# 5 

#Q 

Method 

Objeetive 

Estimated 

Eower Bound 

Run Time (mins) 

7 

1 

Regular 

104 

103 

<1 

Two-step 

104 

103 

<1 

7 

10 

Regular 

104 

103 

<1 

Two-step 

103 

103 

<1 

10 

10 

Regular 

102 

102 

<1 

Two-step 

102 

102 

<1 

15 

10 

Regular 

103 

102 

<1 

Two-step 

103 

102 

<1 

15 

15 

Regular 

103 

102 

<1 

Two-step 

103 

102 

<1 

20 

10 

Regular 

103 

102 

<1 

Two-step 

103 

102 

<1 


5. Norfolk (NMC) 

Norfolk (NMC) is a mid-size ease with 1,189 items. We run SIOM using three 
methods: full, grouping by 500 and grouping by 1,000. (Note: grouping is done in the 
sequential order that the NIINs appear in the input file.) All the methods give the same 
objeetive funetion value of 212. The full model takes 92 minutes to run whereas the 
grouping methods requires mueh less time (25 and 26 minutes respeetively). (See Table 
9.) 


Table 9. Norfolk (NMC) SIOM Results 


Time Eimit 

Group Size 

Objeetive 

Eower Bound 

Run Time (mins) 

2 hours 

500 

212 

211* 

25 

1,000 

212 

211* 

26 

full 

212 

209 

92 


* Estimated, not neeessarily valid. 


The objeetive funetion values suggested by SIOMsQ are slightly above that of 
SIOM for all (s, Q) sets, but the run times of SIOMsQ are mueh smaller. (See Table 10.) 
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Table 10. Norfolk (NMC) SIOMsQ Results 


# 5 

#Q 

Method 

Objective 

Estimated 

Eower Bound 

Run Time (mins) 

7 

1 

Regular 

216 

215 

<1 

Two-step 

216 

215 

<1 

7 

10 

Regular 

216 

215 

1 

Two-step 

217 

215 

1 

10 

10 

Regular 

215 

215 

<1 

Two-step 

215 

215 

1 

15 

10 

Regular 

215 

215 

2 

Two-step 

215 

215 

<1 

15 

15 

Regular 

216 

215 

1 

Two-step 

215 

215 

3 

20 

10 

Regular 

215 

215 

2 

Two-step 

215 

215 

1 


6. Jacksonville (BP28) 

There are 3,256 NIINs in the inventory of Jacksonville (BP28). SIOM is run as 
full and with grouping method restricted to 2 hours and 5 hours for each method. In none 
of these cases is the optimality gap reduced below 1%. The model with time limit of 5 
hours gives the best objective solution (37). (See Table 11.) 


Table 11. Jacksonville (BP28) SIOM Results 


Time Eimit 

Group Size 

Objective 

Eower Bound 

Run Time (mins) 

2 hours 

500 

39 

30* 

122 

1,000 

46 

26* 

122 

full 

121 

25 

122 

5 hours 

500 

37 

32* 

302 

1,000 

37 

27* 

302 

full 

38 

25 

302 


* Estimated, not necessarily valid. 


When we run SIOMsQ, we achieve very reasonable objective function values in 
much shorter time comparing to SIOM. The run time of the (20, 10) case is only 3 
minutes and the achieved objective function value is better than SIOM. 
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As expected, as we increase the number of candidate ( 5 , Q) pairs, we observe 
lower objective function values. Both the regular and two-step methods do not make any 
difference in this case. (See Table 12.) 


Table 12. Jacksonville (BP28) SIOMsQ Results 


# 5 

#Q 

Method 

Objective 

Estimated 

Eower Bound 

Run Time (mins) 

7 

1 

Regular 

41 

41 

1 

Two-step 

41 

41 

1 

7 

10 

Regular 

40 

40 

1 

Two-step 

40 

40 

1 

10 

10 

Regular 

36 

36 

2 

Two-step 

36 

36 

2 

15 

10 

Regular 

34 

34 

2 

Two-step 

34 

34 

2 

15 

15 

Regular 

34 

34 

3 

Two-step 

34 

34 

3 

20 

10 

Regular 

33 

33 

3 

Two-step 

33 

33 

3 


7. Kings Bay (BP28) 

Kings Bay (BP28) has 4,513 NIINs in its inventory. A full and grouped SIOM 
give almost the same objective function value but the running time of grouping heuristics 
are much shorter than the full SIOM. (See Table 13.) 


Table 13. Kings Bay (BP28) SIOM Results 


Time Eimit 

Group Size 

Objective 

Eower Bound 

Run Time (mins) 

2 hours 

500 

132 

132* 

34 

1,000 

133 

132* 

34 

full 

132 

131 

121 


* Estimated, not necessarily valid. 


While running SIOMsQ, and for any (5, Q) pairs generated, the objective function 
value does not change. It is barely greater than the result of SIOM, but run times are 
shorter. 
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For (15, 15) and (20, 10) cases, executing regular SIOMsQ takes more time than 
the two-step SIOMsQ. The main reason for this is that the calculation of an initial 
solution takes very long without using the restricted SIOMsQ'. (See Table 14.) 


Table 14. Kings Bay (BP28) SIOMsQ Results 


# 5 

#Q 

Method 

Objective 

Estimated 

Eower Bound 

Run Time (mins) 

7 

1 

Regular 

134 

134 

3 

Two-step 

134 

134 

3 

7 

10 

Regular 

134 

134 

4 

Two-step 

134 

134 

5 

10 

10 

Regular 

134 

134 

7 

Two-step 

134 

134 

4 

15 

10 

Regular 

134 

134 

4 

Two-step 

134 

134 

4 

15 

15 

Regular 

134 

134 

29 

Two-step 

134 

134 

6 

20 

10 

Regular 

134 

134 

10 

Two-step 

134 

134 

6 


8. Norfolk (BP28) 

Norfolk (BP28) is a large test case with 7,556 items. When full SIOM is run with 
2 and 5 hours limit, the suggested solutions are not acceptable since the differences 
between achieved objective function value and the lower bound are too big. When a 
grouping heuristic is used, the results are much better for both 2-hour and 5-hour cases, 
but still not good enough to determine the answer to the problem. (See Table 15.) 


Table 15. Norfolk (BP28) SIOM Results 


Time Eimit 

Group Size 

Objective 

Eower Bound 

Run Time (mins) 

2 hours 

500 

354 

64* 

127 

1,000 

348 

60* 

127 

full 

17,533 

59 

128 

5 hours 

500 

170 

70* 

307 

1,000 

194 

61* 

307 

full 

13,423 

59 

308 


* Estimated, not necessarily valid. 
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The solutions given by SIOMsQ are better than those by SIOM, since the 
objective function values are close to the provable lower bound of SIOM, and run times 
are acceptable even for larger number of candidates (10 - 20 minutes). The results 
suggested by the (15, 10), (15, 15) and (20, 10) cases are better than those suggested by 
other cases. (See Table 16.) 


Table 16. Norfolk (BP28) SIOMsQ Results 


# 5 

#Q 

Method 

Objective 

Estimated 

Lower Bound 

Run Time (mins) 

7 

1 

Regular 

138 

138 

6 

Two-step 

138 

138 

6 

7 

10 

Regular 

130 

130 

8 

Two-step 

130 

130 

9 

10 

10 

Regular 

99 

99 

12 

Two-step 

99 

99 

10 

15 

10 

Regular 

84 

84 

12 

Two-step 

84 

84 

12 

15 

15 

Regular 

84 

83 

20 

Two-step 

84 

83 

15 

20 

10 

Regular 

80 

80 

14 

Two-step 

80 

80 

13 


9. Bangor (BP28) 

Bangor (BP28) is the second largest test case with an inventory that consists of 
8,141 NIINs. When full SIOM is run, it cannot reduce the optimality relative tolerance to 
1% in the given time limits (2 or 5 hours). The result of full SIOM is not even close to the 
optimal solution, but using the grouping heuristic gives reasonable approximations (an 
objective function value of 287 in the best case, for a lower bound of 251) especially 
when we allow it to run for 5 hours. (See Table 17.) 
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Table 17. Bangor (BP28) SIOM Results 


Time Limit 

Group Size 

Objective 

Lower Bound 

Run Time (mins) 

2 hours 

500 

361 

272* 

129 

1,000 

390 

258* 

129 

full 

10,219 

251 

129 

5 hours 

500 

292 

276* 

309 

1,000 

287 

263* 

309 

full 

6,104 

251 

309 


* Estimated, not neeessarily valid. 


The objeetive funetion values suggested by SIOMsQ are aeeeptable (between 269 
and 271), improving SIOM, and running times are shorter. The number of (s, Q) pairs or 
the two-step method does not make a signifieant differenee in solution quality for this 
ease. (See Table 18.) 


Table 18. Bangor (BP28) SIOMsQ Results 


# 5 

#Q 

Method 

Objective 

Estimated 

Lower Bound 

Run Time (mins) 

7 

7 

Regular 

270 

270 

8 

Two-step 

270 

270 

8 

7 

10 

Regular 

271 

270 

12 

Two-step 

271 

270 

14 

10 

10 

Regular 

269 

269 

36 

Two-step 

269 

269 

28 

15 

10 

Regular 

270 

269 

15 

Two-step 

270 

269 

14 

15 

15 

Regular 

271 

269 

32 

Two-step 

271 

269 

32 

20 

10 

Regular 

271 

269 

17 

Two-step 

271 

269 

17 


10. Yokosuka (BP28) 

There are 11,798 items in the Yokosuka (BP28) test ease. Due to its large seale, 
the SIOM solution quality is not good enough sinee the optimality relative toleranee 
eannot be redueed below 1% even if grouping heuristie is used. The best solution is 
achieved for the 5 hour - 500 group size run with an objective function value of 4,211 and 
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a (provable) lower bound of 3,321 (i.e., a potential solution gap of at most 26.8%). (See 
Table 19.) 


Table 19. Yokosuka (BP28) SIOM Results 


Time Limit 

Group Size 

Objective 

Lower Bound 

Run Time (mins) 

2 hours 

500 

5,093 

3,393* 

138 

1,000 

15,335 

3,391* 

137 

full 

45,820 

3,310 

137 

5 hours 

500 

4,211 

3,433* 

318 

1,000 

8,191 

3,399* 

318 

full 

32,401 

3,321 

322 


* Estimated, not necessarily valid. 


When we run SIOMsQ, we obtain acceptable solutions in a reasonable amount of 
time. Specifically, we obtain an objective function value under 3,700 in all runs, which 
means much smaller solution gaps (always under 11.4%). This calculation uses the 
provable lower bound of SIOM, given SIOMsQ’s bound is only an estimate. (See Table 
20 .) 

There are 2,041,029 integer variables and 82,589 single equations in SIOMsQ run 
for the (15, 15) case in Yokosuka (BP28). 


Table 20. Yokosuka (BP28) SIOMsQ Results 


# s 

#Q 

Method 

Objective 

Lower Bound 

Run Time (mins) 

1 

1 

Regular 

3,686 

3,661 

30 

Two-step 

3,686 

3,661 

30 

1 

10 

Regular 

3,661 

3,661 

118 

Two-step 

3,691 

3,661 

58 

10 

10 

Regular 

3,664 

3,660 

81 

Two-step 

3,665* 

3,660* 

97* 

15 

10 

Regular 

3,659 

3,659 

93 

Two-step 

3,659 

3,659 

91 

15 

15 

Regular 

3,661 

3,659 

171 

Two-step 

3,660 

3,659 

173 

20 

10 

Regular 

3,684 

3,658 

135 

Two-step 

3,663 

3,658 

140 


* For hi = 30 (instead of /?' = 50 used in all other cases) 
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C. SUMMARY OF RESULTS 

The results of the all test eases are summarized in Table 21. We present the best 
solutions aehieved by SIOM along with its lower bounds and running times for each test 
case. We also show the solutions, candidate pair sizes, execution times and the solution 
methods of best SIOMsQ for every case. 

For all of the test cases, the solutions that SIOMsQ achieves are better than or 
very close to that of SIOM, but the running times of SIOMsQ are substantially shorter. 
As the size of the test cases increases, the best SIOMsQ solutions are achieved with the 
large candidate sizes such as (10, 10), (15, 15) or (20, 10). The two-step approach 
(including the time spent in solving the first step model, SIOMsQ') is faster especially in 
the larger test cases. 
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Table 21. Results Summary of Test Cases 


Test Case 

Best SIOM 
Solution 

Provable SIOM 
Lower Bound 

SIOM Time 
(mins) 

Best SIOMsQ 

#Q) 

Best SIOMsQ 
Solution 

Best SIOMsQ 
Time (mins) 

Best SIOMsQ 
Method 

Yokosuka 

(NMC) 

0 

0 

<I 

Multiple 

0 

<I 

Both 

Key West 
(BP28) 

1.40 

1.39 

<I 

(15, 15) 

1.38* 

<I 

Both 

Kings Bay 
(NMC) 

30.7 

30.7 

<I 

Multiple 

30.7 

<I 

Both 

San Diego 
(NMC) 

103 

102 

19 

(10, 10) 

102 

<I 

Both 

Norfolk 

(NMC) 

212 

209 

92 

Multiple 

215 

<I 

Both 

Jaeksonville 

(BP28) 

37 

25 

302 

(20, 10) 

33 

3 

Both 

Kings Bay 
(BP28) 

132 

I3I 

I2I 

Multiple 

134 

3 

Both 

Norfolk 

(BP28) 

170 

59 

307 

(20, 10) 

80 

13 

Two-step 

Bangor 

(BP28) 

287 

251 

309 

(10, 10) 

269 

28 

Two-step 

Yokosuka 

(BP28) 

4,211 

3,321 

318 

(15, 10) 

3,659 

91 

Two-step 


* Due to a better eoverage in Q, the optimal solution by SIOMsQ in the (15, 15) ease ean potentially be better than that of SIOM. 
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IV. CONCLUSIONS AND FUTURE RESEARCH 


A. CONCLUSIONS 

In this thesis, we have developed SIOMsQ, a new approaeh to solving SIOM (a 
mathematieal program developed at the Naval Postgraduate Sehool to guide NAVSUP’s 
inventory levels of site demand-based items). 

We have sueeessfully implemented SIOMsQ, whieh restriets SIOM by redueing 
the number of eandidate order-points in exehange for more effieient solutions. The 
approaeh ereates an edueated speeifieation of a list of eandidate (y Q) pairs and solves 
the optimization problem using only those eandidates as deeision variables. SIOMsQ may 
still be diffieult to solve in some eases, so we further restriet it using a two-step approaeh 
in order to find an initial solution with whieh to speed SIOMsQ up. 

We run 10 test eases provided by NAVSUP with various eombinations of 
eandidate pair sizes. In eight of ten eases, we obtain at least as good solutions as SIOM in 
the same or substantially less time. In the two other eases, the differenees are minimal, 
and eomputational times are still favorable. 

As a eonelusion of this researeh, we reeommend that running SIOMsQ using 
(#5, #Q) values of (20, 10) or (7, 7), and a two-step approaeh are the preferred options. 
(7, 7) oases give good-enough results faster than (20, 10) oases, but the solutions 
suggested by (20, 10) oases are more preoise. Therefore, there is a tradeoff between 
solution quality and solving time, although in most oases, differenees are also small. We 
reeommend using SIOMsQ versus the ourrent SIOM. 

B. FUTURE RESEARCH 

As we have suooessfully implemented SIOMsQ, further studies oan foous on: 

• Having SIOMsQ fully tested by NAVSUP, 

• Integrating SIOMsQ as an option for NAVSUP to use from the graphioal 
interfaoe of SIOM, 
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• Analyzing the optimal configuration of (#s, #Q) pairs for each individual 
test case, or even developing a methodology that allows SIOMsQ to 
dynamically accommodate different candidate pair sizes. 

• Speeding up pre-processing of generating (5, Q) candidate pairs, the 
probability distributions of demand based on Q, and derived data based on 
5 and Q (all of which takes up to 50 minutes in some cases). 
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APPENDIX. SIOM FORMULATION 


This appendix describes the formulation of SIOM as given in Salmeron and 

Craparo (2016). 

A. INDICES AND INDEX SETS 

* , item (i.e., NUN), for i e /. 

^ , Group. (Note; I here could also represent a combination of group and other 
criteria, such as group A, B, C, D, for the purpose of establishing different 
target fill rates.) 

q, index for candidate order quantities, for q e {1,...,^.} (see q. below). 

” , demand-level index, for n&N (e.g., N = {1,2,...,10} represents ten levels 

of demand). This is actually item-dependent, and order size dependent, see 
n‘‘ parameter below. 

^ , penalty segment, for msM (e.g., M = {1,2,...,5} represents five levels 
of penalties for deviations, with respect to desired fill rate levels). 

B. INPUT DATA AND PARAMETERS 

t ., lead time length for item i [quarters/lead time]. 

X;, , expected demand and estimated deviation, respectively, for item i during 

the lead time [items/lead time], [items/lead time]. 

q ^, number of candidate order quantities for item i . This number depends on 
X,., ,and 5,.. 

gf , q-ih candidate order quantity for item i [units/order]. 
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cycles during a lead time for item i and q-th eandidate order quantity. 


c‘‘ 


Caleulatedas cf =max{l,^} [orders per lead time]. 


, ealeulated as: one if x,. +1 > , and zero otherwise. 


A[, ealeulated as: A[ = max{A[^}. 

q\q<qi 


expeeted demand and estimated deviation adjusted for eyele time, for 
q-th eandidate order quantity for item i. Caleulated as: = x,. / c. and 

0^=0'/c:^. 

fi, desired (target) fdl rate for eaeh item in group / [fraetion]. 


Wi , weight for meeting required fdl rate eaeh item in group / [weight units]. 


, penalty for deviating from required fdl rate for items in group I within 
penalty segment m . Caleulated as w,^=m‘' w,, where e is a user input 
(penalty exponent, default e = 1). 

fii„, maximum deviation allowed for fill rate for item/ within penalty segment 


m . 


Caleulated as = fi^ 


m 

j^M 


where / e . 


OL., operating level for NUN / [months]. 
SLi, safety level for NIIN / [months]. 


5, , allowanee for item / [demand units]. 


A-, one if allowanees are aetivated, and zero otherwise, 
q, eost per unit at MSL [$ / unit]. 
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b. 


overall budget [$]. Caleulated as: 


i 


c, 


max 10, '^OL,+SL;)xf — 


^ L (OL,+5L,)xf 1 

max 1 0, ^-Li_!_ + X, 

3 


, if allowanees are aetivated (A- = 1), or 

otherwise (A- = 0), 


where =average quarterly demand for NUN / [units/quarter], derived as 


= X; /1,.. 


p2,d[^ , n-th level of probability and demand, respeetively, adjusted for the eyele 
time for the q-\h eandidate order quantity for item i [unitless], [demand 
units]. Note: If probability distributions are estimated based on mean and 

varianee, use x[‘‘ and to generate these. 


n‘f, number of demand levels, adjusted for the eyele time for the q-ih 
eandidate order quantity for item i. 

derived data: = max^^^, {dlP^ ; [demand units]. 

r , maximum number of total expeeted orders per month [orders]. Note: Can 
be fraetional. 


lower and upper bounds on reorder point for item i [items], [items]. 
A >-1 . 

, shelf life for item i (to establish its MSL upper bound) [months]. 

S , maximum months of supply for any item (to establish its MSL upper 
bound) [months]. 

initial reorder point used to enforee persistenee for item i [items]. 

Sf , penalty for deviation from initial reorder point for item i [weight units]. 


35 



y’’, overall persistence level [weight units]. 

, penalty for item’s i deviation above maximum months of supply [weight 
units]. 

§ 

y , overall weight for MSL deviations above maximum months of supply 
[weight units]. 

, large number greater (in magnitude) than any possible ‘negative fill rate’ 
value given by the original fill rate formula for the <y-th candidate order 
quantity for item i, if = 1. For example, we may use Mf = 1,000. 
However, we can estimate a better value for numerics: 
= 2- x'VQ^ l/e« 

i i i 

M :, calculated as M = max {Mf}. 

C. DECISION VARIABLES 

5;, reorder point for item i [items]. 

2,., order quantity for item i [items]. 

Ff, one if the <y-th candidate order quantity for item i is selected, and zero 
otherwise. 

, deviation (down and up, respectively) with respect to initial reorder point 
for item i [items]. 

5/, deviation of item’s i MSL above the average demand during the 
maximum months of supply for the item [months]. 

fi ,fi , fill rate and negative component of fill rate (if any), respectively, for item 
i [fraction]. (The negative component is only applicable if c. > 1.) 
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, deficit in fill rate (with respect to target) for item i in penalty segment m 


[fraction], 

, ancillary variable for stockouts for item i , given its q-ih candidate order 
quantity is chosen and the demand level 

n<n1: = max {<• - (x, - (c* -1)6?), o}. 

ancillary binary variables to signal the fill rate sign, if cf >\. 


D. FORMULATION 


min 

s,Q,r,s^,s“,f,f,z 




Ij m J im I 


(B) 


subject to: 



e, = Z ST? Vi 


(14) 

Zrf = i Vi 

9|?£'7, 


(15) 

n\n<nf 

\fi,q\q<q^ 

(16) 

n\n<nf 

VI 

> 

(17) 

Vf|A[=l 


(18) 

f->fr/M, Vi|Af=l 


(19) 

f,+fr=l Vi|Af=l 


(20) 

4’®° ^ - (^i - (4 - l)2f ) yi,n,q\n<n‘l,q<q. 


(21) 
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(22) 




1 

+ 


(23) 

i ?!?<?, 

■r»<r 

(24) 

rslA^Si,+e,< x,"!- Vi 

(25) 

^(i, +a)<s+s* Vi 

it/ 

(26) 

^i = + sl-s 

T Vi 

(27) 

^.<5. <5^ 

Vi (optional) 

(28) 

Q. > 0 and integer V/ 

(29) 

5; > -1 and integer Vi 

(30) 

He {0,1} 

Vi,^k<^, 

(31) 

~^i 

r^i 

+ 

IV 

O 

Vi 

(32) 

?,so>o 

^in 

yi,q,n\n<nf ,q<q^ 

(33) 

o< /."</- 

^ im J im 

\fi,m 

(34) 

01 

1 

Vi|A;=l 

(35) 


Vi 

(36) 
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